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Abstract. Production and elliptic flow of heavy quarks are investigated in nucleus-nucleus 
collisions at RHIC and LHC within the Boltzmann Approach of MultiParton Scatterings 
(BAMPS). The initial heavy quark yield is estimated with the leading order mini-jet model 
O l' and compared to PYTHIA for several parton distribution functions. Secondary production of 

D , heavy quarks in the quark-gluon plasma is examined within the full 3+1 dimensional BAMPS 

43 i simulation of heavy ion collisions. At RHIC this yield is negligible, but for LHC charm 

production in the QGP plays a significant role. In addition, we study the elliptic flow of charm 
quarks at RHIC with the result that leading order processes are not sufficient to describe the 
experimental data. 
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■^j- ■ 1. Introduction 

Charm and bottom quarks are an interesting probe for the hot and dense medium, which is 
produced in heavy ion collisions at RHIC [JJ [5J El H] . Due to their large mass heavy quarks are 
produced at a very early stage of the collision, revealing unique insight into the initial properties 
of the quark-gluon plasma (QGP). This large mass also determines the relevant scale of their 
^ ■ production processes, which should, therefore, be describable in the framework of perturbative 

QCD (pQCD). 

The production of heavy quarks can be divided in three stages: They are created (i) in 
hard scatterings during initial nucleon-nucleon collisions, (ii) in the QGP and (hi) during the 
hadronic phase. We will neglect the latter one since the energy density in the hadronic gas is low 
and heavy flavor production becomes very unlikely. At RHIC and LHC primary heavy quark 
production in hard scatterings plays the most important role. At RHIC secondary production 
in the QGP is nearly negligible, but at LHC it contributes significantly to the total heavy flavor 
yield. 

The experimentally observed elliptic flow [5] and energy loss [61 [5] of heavy quarks indicate 
that they interact strongly with the medium and thermalize relatively fast, although the "dead 
cone effect" [8] predicts a delay in their thermalization compared to light quarks. The 
investigation of this puzzle is still in progress [9j [TOj, [Til 021 E2J G3] ■ 

The present article is organized as follows: In the next section we describe our model, the 
parton cascade BAMPS. Thereafter, we discuss primary heavy quark production in hard parton 
scatterings during initial nucleon-nucleon collisions within the leading order mini-jet model 
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and PYTHIA. In ISec, 41 secondary heavy quark production in the QGP at RHIC and LHC 
is investigated with BAMPS for different initial conditions. Finally, we address elliptic flow of 
charm quarks at RHIC in ISec. 5} followed by brief conclusions. 



2. Parton cascade BAMPS 

The full 3 + 1 space-time evolution of the QGP is studied within the partonic transport model 
Boltzmann Approach of MultiParton Scatterings (BAMPS) [15(. I16j . which solves the Boltzmann 
equation, 

(s + !s)«™ t >= < ^ ,+ *" ,+ -' (1) 

dynamically for on-shell partons with a stochastic transport algorithm and pQCD interactions. 
Ci denotes the collision integrals of the 2 — > 2 and also 2 3 interactions and /j(r, Pi,t) the 
one particle distribution function of species i = g, Q, Q (Q = c,b), since light quarks are not 
included yet. The following processes are implemented in BAMPS: 

9 + 9^ 9+9 
9+9^9+9+9 
9+9+9^9+9 
g + g^Q + Q 
Q + Q -> 9 + 9 
9 + Q -> 9 + Q 

g + Q^g + Q (2) 

For further details on the purely gluonic processes we refer to [151 US] an d to [T7] for the 
processes, in which heavy quarks are involved. 



3. Primary heavy quark production in initial nucleon-nucleon scatterings 

For the estimation of the primary heavy quark yield during initial nucleon-nucleon scatterings 
we use the following two approaches. 



3. 1 . Leading order heavy quark production within the mini- jet model 
In leading order (LO) heavy quarks are produced in the two processes 

g + g^Q + Q 

q + q^Q + Q . (3) 



The differential cross section for heavy quark production in the collision of hadron A and 
B is essentially given by the product of the parton distribution functions (PDF) /j and the 
differential partonic cross sections for both processes, [18] 



dd 
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(4) 



y being the rapidity and pt the transversal momentum of the heavy quark and anti-quark in 
the center of mass frame. The PDFs are dependent on the factorization scale and the partonic 
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PDF 


Reference 


Charm pairs 


CTEQ51 (LO) 


[23j 


8.9 


CTEQ61 (LO) 




M 


9.2 


CTEQ6m {MS) 


[24j 


13.6 


MRST2007LOmod 


[25j 


9.2 


HERAPDF01 


[26j 


12.3 


GRV98 (LO) 


[27j 


3.0 



Table 1. Number of charm pairs produced in primary hard scatterings in central Au+Au 
collisions at RHIC for some parton distribution functions by sampling nucleon-nucleon collisions 
with PYTHIA and scaling to Au+Au collisions. 



PDF 


Reference 


Heavy q 
Charm 


uark pairs 
Bottom 


CTEQ61 (LO) 


[24| 


62 


7.2 


CTEQ6m {MS) 


[24| 


66 


6.9 


MRST2007LOmod 


[25] 


67 


8.9 



Table 2. As lTab. II but for central Pb+Pb collisions at LHC. 



cross sections on the renormalization scale. Therefore, heavy quark production in general is 
sensitive to both scales. Other uncertainties come from the heavy quark mass and the choice 
of PDFs. Depending on these parameters the total LO heavy quark production cross section in 
nucleon-nucleono collisions changes by a factor of about 2 at RHIC and even of about 4 at LHC 
energy [17] . 

3.2. Heavy quark production with PYTHIA 

In addition, we use PYTHIA [19] . the LO event generator for nucleon-nucleon collisions, to 
estimate the number of initially produced heavy quarks in heavy ion collisions. For that we 
scale heavy quarks from nucleon-nucleon collisions with roughly the number of binary collisions- 
Details on how this is done can be found in [17J. Tables [1] and [2] list the number of initially 
produced charm quarks in central Au+Au collisions at RHIC and charm and bottom quarks 
in central Pb+Pb collisions at LHC, respectively, according to PYTHIA for various parton 
distribution functions. As a note, there are also large errors due to uncertainties in mass, 
shadowing, factorization and renormalization scale like in the LO pQCD calculation, although 
these are not reflected in the tables. 



Fig. 1 compares the rapidity distribution of the charm production cross section in one nucleon- 
nucleon collision simulated with PYTHIA together with the LO pQCD calculations and the 
experimental data points |28} 129], The LO calculations are far below both experimental values, 
which indicates that next-to-leading order processes [31] have to be considered. In contrast, 
PYTHIA lies much closer at least to the PHENIX data point. Therefore, we decided to use 
PYTHIA for sampling our initial heavy quark distributions. For the PDFs we chose CTEQ61, 
because it is designed for LO event generators such as PYTHIA 



Actually, this value is lowered to 1000 at RHIC [20] and 1500 at LHC [211 122| in order to take shadowing into 
account . 
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Figure 1. Charm production cross section dcr^ /dy as a function of rapidity y in a nucleon- 
nucleon collision at RHIC energy simulated with PYTHIA and pQCD, respectively, for the PDFs 
CTEQ61 and CTEQ6m together with experimental data |28[I29|. The pQCD calculation is done 

in LO with l i F = m = Jp^ + M?, N f = 3, M c = 1.5 GeV, A QC d = 346 MeV [30] and K = 1. 



4. Secondary heavy quark production in the QGP 

4-1- Charm production at RHIC 

We study charm production during the QGP phase in central Au+Au collisions at RHIC within 
the framework of BAMPS. The initial charm distribution is obtained from PYTHIA. Initial 
gluons are sampled according to three different models: a color glass condensate (CGC) inspired 
model [32^ I33j . the mini-jet model \34\ 135] . and PYTHIA. The detailed prescription how we 
scale partons from PYTHIA to heavy ion collisions can be found in |17| . For the partonic cross 
sections in BAMPS of the processes from (2) we employ a constant coupling of a s = 0.3. 



Fig. 2 shows the charm yield in the QGP as a function of time. For all three initial distribution 



models the charm production during the QGP is very small compared to the initial yield. With 
PYTHIA initial conditions only 3 % of the total final charm quarks are produced in the QGP. 
For CGC or mini-jet initial distributions for the gluons this fraction is about 6%. The higher 
yield in the latter models is due to the higher initial energy density [17]. The main increase in 
the charm number takes place at an early stage of the evolution, where the QGP is not fully 
equilibrated yet. 

Changing the charm mass to a smaller value of M c = 1.3 GeV or employing a K factor of 2 
to the charm production cross section increase these results by a factor of 2, respectively. As a 
result, using mini-jet initial conditions, M c = 1.3 GeV and also K = 2 gives the maximum value 
in our model of 3.4 charm pairs produced in the QGP (cf. |Fig. 3 ), which is still just 27% of the 
total charm yield. 

Our results indicate that the charm production in the QGP at RHIC is nearly negligible. That 
is also underlined by experiment, where a charm scaling with the number of binary collisions 
was found [361 [281 [37] . Charm quarks are, therefore, produced in initial hard scatterings during 
nucleon-nucleon collisions. 

At first sight this small charm production is a bit surprising since the charm quark fugacity 
is always smaller than 1 in the QGP [17] , However, the reason for that lies in the huge timescale 
of chemical equilibration of charm quarks in the medium produced at RHIC, which is about a 
factor of 100 larger than the lifetime of the QGP [17J. 
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Figure 2. Number of charm quark pairs produced in a central Au+Au collision at RHIC 
according to BAMPS. The initial parton distribution for gluons is obtained with PYTHIA, CGC, 
and the mini-jet model. For better comparison the initial charm distribution from PYTHIA is 
used for all models. Please note the small range of the y axis. 
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Figure 3. As Fig. 2 varying initial gluon conditions, charm mass M and K factor for gg 
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in addition, the experimental value for the number of final charm pairs is plotted [36] 



4-2. Charm and Bottom production at LHC 

At LHC the picture changes compared to RHIC: Due to the larger energy density, charm 
production in the QGP cannot be neglected anymore. Fig. 4| depicts our predictions for the 



charm number evolution at LHC. Again, we use PYTHIA for the initial charm yield and 
PYTHIA, CGC, and mini-jet initial conditions for the gluons. With PYTHIA initial conditions 
15 % of the total charm quarks are produced in the QGP. With gluons from CGC this fraction 
even rises and for mini-jet initial gluon conditions the charm yield during the QGP phase is 
nearly equal to the initial charm number. Lowering the charm mass or introducing a K factor 
increase these values again. 

Just like at RHIC, the charm quark fugacity at LHC is far below 1 [T7J. However, because of 
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Figure 4. Number of charm pairs in a central Pb+Pb collision at LHC simulated with BAMPS 
for PYTHIA, CGC, and mini-jet initial conditions for gluons. Charm quarks are sampled with 
PYTHIA for all three models. 



the much smaller equilibration time of charm quarks in such a hot medium as at LHC, which is 
of the same order of magnitude as the QGP lifetime [T7j, the charm production is much stronger. 

As shown in [IT] the bottom production during the QGP phase at LHC and, of course, also 
at RHIC is completely negligible due to the larger mass. Therefore, bottom quarks are just 
produced in initial hard scatterings, which renders them as a promising probe of the early stage 
of the collision due to the precisely determined production time. 




5. Elliptic flow of charm quarks at RHIC 

The elliptic flow of charm quarks is defined by the average 

(5) 

over all charm quarks at mid-rapidity and is an indicator how strong charm quarks are coupled 
to the bulk medium. In the present article we study only charm scattering with gluons in leading 
order, g + c — > g + c. Without considering higher order corrections or employing a K factor for 
the cross section the interaction between charm quarks and the gluonic medium is rather small. 
Therefore, although a strong gluonic v% is observed in BAMPS [38, 39, 40J, which is in agreement 



with experimental data, the of charm quarks with K = 1 vanishes nearly, as shown in Fig. 5 
For a substantial charm flow one must multiply the cross section of gc — > gc by at least a factor 
of 10. In order to come in the region of the experimentally observed flow of heavy flavor [SJ HI] 
even an unreasonable K = 50 is necessary, which shows that the charm scattering cannot be 
described just by the LO process. Of course, one has to be careful with comparing the flow 
from charm quarks to the experimentally observed flow of heavy flavor electrons, which stem 
from the decay of D or B mesons. However, in [32] it is found that the effect of hadronization 
and decay on the flow is negligible. In contrast, in the coalescence model the flow of D mesons 
is considerably larger than that of charm quarks [33]. We also made some calculations with 
bottom quarks and saw that our here presented results are only slightly modified and that only 
for high pr, if one considers both charm and bottom. 
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Figure 5. Elliptic flow of charm quarks with pseudo-rapidity \r]\ < 0.35 at time t = 5fm/c for 
RHIC collisions with impact parameter of b = 8.2 fm. The cross section of gc —> gc is multiplied 
with different K factors. 

In order to improve our calculations we will investigate the effect of a running coupling and 
modify the evaluation of the screening mass in the gc — > gc cross section |444 [131 fl^l I45j . In 
addition, we want to study the impact of taking higher order corrections like g + c — > g + c + g 
into account, for which BAMPS is an ideal framework since 2 3 interactions are already 
included for gluons |15j . With these improvements we will also extend recent studies on the 
nuclear modification factor [4<j\ |4"7] to the heavy flavor sector. 

6. Conclusions 

We have studied heavy quark production in heavy ion collisions at RHIC and LHC. The initial 
yield was estimated with LO pQCD processes, which did not agree with the data, and PYTHIA. 
The latter was then used to sample the initial heavy quark distribution for the QGP phase, 
which was simulated with BAMPS. In that framework we investigated the dependence of heavy 
quark production on different initial gluon distributions, namely from PYTHIA, CGC and the 
mini-jet model. 

At RHIC the production of charm and, of course, bottom quarks in the QGP is highly 
suppressed due to the huge chemical equilibration time scale in the medium created at RHIC. 
For LHC we predict a significant charm production of 10 to 60 charm pairs in the QGP phase, 
depending on the initial conditions of the medium. This is due to the high initial energy density 
and a relatively small chemical equilibration time scale, which is of the same order as the QGP 
lifetime. In contrast, bottom production in the QGP at LHC is negligible due to the high bottom 
mass. 

The experimentally observed elliptic flow of heavy quarks cannot be described with LO pQCD 
cross sections. Therefore, we want to take higher orders into account and study the effect of a 
running coupling and improved Debye screening. 
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